Connection between the "Strutinsky level density" and the semiclassical level density 
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We establish an analytical link between the level density obtained by means of the Strutinsky aver- 
aging method, and the semiclassical level density. This link occurs only in the so-called "asymptotic 
limit" . 

It turns out that the Strutinsky method amounts to an approximation to the semiclassical method. 
This approximation contains an unavoidable remainder which constitutes an intrinsic noise in com- 
parison to the semiclassical method. Thus, the "old" problem of the dependency of the Strutinsky 
procedure on the two free smoothing parameters of the averaging is intimately connected to this 
noise. 

On the other hand, we demonstrate that the noise of the method is small in the average density 
of states and in the average energy, whereas it might be non-negligible in the shell correction itself. 
In order to improve this method, we give a "rule" which consists simply of minimizing the relative 
error made on the average energy. 



PACS numbers: 21.10Dr, 21.10.Ma, 21.60.-n 

I. INTRODUCTION 

The inclusion of the Strutinsky's shell correction 
3 in the liquid- drop model [3, namely the so-called 
macroscopic-microscopic method, has allowed consider- 
able improvements in the predictions of the nuclear 
masses [5|, and in the calculations of the fission barri- 
ers, as well 010- Nowadays, despite the progress of the 
more basic microscopic models (such as self-consistent 
models), it remains often used. 

This method consists essentially of combining the 
liquid-drop model (macroscopic model) where the bind- 
ing energy varies slowly as a function of the number of nu- 
cleons N and Z, with a shell correction varying abruptly 
with N and Z . The latter is due to the non-uniformity 
of the shell structure of the energy levels. It is extracted 
from a single-particle Hamiltonian (microscopic model) 
according to an original idea of Strutinsky. 

The Strutinsky method is mainly based on a particular 
smoothing procedure of the density of states. 

Although this method is very efficient, it contains two 
weak points which are: 

1) The dependence of the results on the two well known 
inherent parameters, i.e. the width 7 and the order M 
of the smoothing 

2) The difficulties of the treatment of the continuum 
encountered with realistic mean potentials 

The purpose of the present work is summarized in the 
following points: 

1) The Strutinsky method can be derived rigorously 



from the point of view of the least-squares approximation 
of the level density. The equivalence between this point 
of view and the well-known standard averaging appears 
trivial. 

2) In this work, it is proved analytically that the aver- 
aged level density obtained by the Strutinsky method is 
simply an approximation to the semiclassical level den- 
sity. 

In this respect, the semiclassical level density can be con- 
sidered as the "true" (i.e. exact) smooth density. 

3) It is mainly shown that in comparison with the semi- 
classical method, the Strutinsky method is characterized 
by a remainder which contains all the dependence on the 
two smoothing (free) parameters, and which is thence, 
the source of the "noise" of the averaging procedure. 

4) Concerning the smooth density of states and the 
smooth energy, it is demonstrated that the Strutinsky 
method is reliable. However, the shell correction itself 
must be treated with care because it is very sensitive to 
the choice of the two free averaging parameters. In this 
context, in order to improve the method, we propose the 
"rule of the relative remainder" (RRR rule). 

5) At last, it is explained why the Strutinsky method 
fails near the zero-energy (top of the well) for finite mean 
potentials 



II. THE STRUTINSKY AVERAGING 

A. Bases and phenomenology of the Strutinsky's 
method 
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In spite of the complexity of the nuclear forces, it ap- 
pears that most of the binding energy of the nuclei is 
well described by the so-called liquid-drop model. This 
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simple phenomenological approach is of a classical type. 
This means that the quantum effects, or a more precisely, 
the shell effects, are ignored by this model. This causes 
systematic discrepancies between the theoretical predic- 
tions and the experimental data 

On the other hand, it is known that such effects are con- 
tained in the shell model, but the latter is unable to re- 
produce correctly the general trends of the binding en- 
ergy. To solve this dilemma, Strutinsky has proposed 
to combine the binding energy of the liquid-drop model 
with a small (but essential) correction deduced from the 
shell model. This can be written as: 



ref. the authors use a method referred to as the Ex- 
tended Thomas Fermi plus Strutinsky integral (ETFSI) 
method. In the latter, the semiclassical quantity i?(p) is 
deduced self-consistently from a microscopic effective in- 
teraction. (Skyrmc III). The shell correction is then added 
to E{p). 

In the following, we will mainly focus on two points: 

• The mathematical aspect of the Strutinsky smooth- 
ing 

• The link between this smoothing and the semiclas- 
sical approximation, and its consequences 



E (Binding Energy) 



E (Liquid Drop Model) 
+5E (Shell Correction) 



The shell correction is calculated from a mathematical 
prescription outlined by Strutinsky 

It is obtained by summing the single particle energies of 
a phenomenological shell model potential and subtracting 
the average (smooth) part of this quantity: 



SE (Shell Correction) = — ei 



As already mentioned, this method is often called the 
macroscopic-microscopic method because it mixes two 
very different models. Such a duality is obviously not free 
from inconsistency. Nevertheless, it is possible to give a 
microscopic basis to this "model" within the Hartree- 
Fock (HE) approximation and some simple assumptions 

a,i. 

This consists essentially of expanding the HF energy 
around its semiclassical approximation obtaining thus 
the so-called " Strutinsky energy theorem" : 



E{p)=Eip)+[Y,e,-J2, 



B. The "exact" level density 

The Strutinsky averaging can be derived by various 
formal approaches [f, 

In this work, we will derive it from the point of view of 
the least-squares approximation. In fact, this point has 
been early suggested by Bunatian and co-authors [l^. 
The least-squares approximation will help us to under- 
stand why the Strutinsky averaging fails near the edge of 
realistic potentials. 

For an entirely discrete spectrum the level density of 
states is defined by: 



(2) 



O2 



(1) 



In the following, this density will be called the "exact 
quantum level density" , or shortly, the " exact level den- 
sity" , because it is a true quantum quantity, as opposed 
to its semiclassical approximation, or as opposed to the 
level density obtained from (70(e) by the Strutinsky aver- 
aging (see below). 

In fact, Eq.(I2]) concerns uniquely infinite potentials 
without continuum. Finite potentials will be treated sep- 
arately at the end of this paper. 



Here p is the HF density matrix, and p its semiclas- 
sical approximation which is a smooth quantity, free of 
shell effects. For this reason p can be assimilated to the 
classical average (i.e. without quantal variations) part of 
p. The sums of single-particle energies Ee^ and Se^ are 
related respectively to p and p. Finally, O2 is a quantity 
of the second order in the operator p — p, and is generally 
negligible (for details see [28l|). 

It is clear from (jT]) that the macroscopic-microscopic 
method described above is a "schematic" interpretation 
of this theorem. The quantity E(j)) can be replaced by 
the energy of the liquid-drop formula, and the shell cor- 
rection is deduced from a phenomenological one body 
Hamiltonian 

It is also to be noted that a complete microscopic ap- 
proach of this theorem remains possible. For example, in 



C. Polynomial approximation to the "exact" level 
density: 

Let (7a/ (e) be a polynomial approximation of order M 
to the "exact" level density. More precisely, we seek 
this approximation in the vicinity of a point A (which 
represents actually the Fermi level) in an effective in- 
terval [—7 -I- A, A -I- 7], by using the Gaussian weight 
exp (— (e — A)^/7^) . For this reason, the cited polyno- 
mial must depend a priori not only on Af, but also on 7 
and A. Therefore, it will be denoted gj^,j ^(e. A). For our 
purpose, it will be convenient to write this polynomial as 
a linear combination of Hermite polynomials Hk- 



M 



9m 



^(e,A) = ^Cfei/fc 



fe=0 



e- A 
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Thus, We must look for the polynomial (7M,7(e, A) 
which minimizes the integral: 



e- A 



/(A,M,7)= / [9o{e)-gM e \ ^ J de 



(3) 

This procedure is a local averaging in the sense of the 
least- squares fit. Minimizing ([3]) with respect to the 
coefBcients Ck , and using the orthogonality property of 
the Hermitc polynomials we find: 



M 



,(e,A)= ^c„,(A,7)i?,n(^) (4) 



9m 



m=0 
oo 



1 1 

Cm(A, 7) = V i/„(w„)- exp [-ul] , 



n=0 

£^^ - A 

7 



As we shall see in the next subsection, Eq.lH]) is not 
the final definition of the averaging which thus appears 
somewhat more " subtle" . 



9m , = Xl^A/ (Un) , 



(5) 



Fm {x) = Pm {x) - exp (- x^) 



1 

7 



M 



m=Q V 

H„,{0) ^ (-l)"'/2.m!/ (m/2)! if m is even, 
H^{0) = if 771 is odd 

In ^ the polynomial Pm constitutes the so-called cur- 
vature correction term. 

It is easy to check that the expression ([5]) obtained 
from a least-squares fitting, can be written as the usual 
folding procedure of the exact level density , [3 ■ 



9,a^) 



9o{e)FM 



7 



dt 



(6) 



which demonstrates the equivalence between the two 
points of view, i.e. between the local least square smooth- 
ing dS]) and the averaging 
Remark: 

It is to be noted that if M increases to infinity and/or 7 
decreases to zero, then 'g^^^^ tends toward go(e)- Obvi- 
ously, in practice these parameters are finite. 



E. Necessary condition in the smoothing procedure 



D. The Strutinsky's averaging as a moving average 

In Eq.jH) the Fermi level A is supposed to be fixed, 
and the polynomial 'g^ ^(e, A) smoothes the exact level 
density only in the vicinity of A. Thus only a part of the 
spectrum is " smoothed", i.e. the part e ~ A. To avoid 
this drawback, it is necessary to consider A as a variable. 
The averaged level density is thus defined as 'g^ ^(e. A) 
making e = A. This amounts to perform a moving aver- 
age (i.e. A is moved with e, "centering" e on A). 
In the following, we shall call up the function gj^j^y{X, A) 
as the Strutinsky's level density, and we will note it sim- 
ply by 5j,,,7(A). 

Since the coefficients Cm in (IH) depend on A, the Strutin- 
sky's quantity -y(A) is, in general, not a polynomial 
in A. However, it is clear that, although the Strutinsky 
level density is not really a polynomial, it behaves locally 
(~ A ± 7) , like a least deviating polynomial approxima- 
tion for the exact level density .90(A) given by ([2|). 

Explicitly, we must replace e with A in ([¥]): 



The least-square smoothing ([3]) or its equivalent ^ 
gives an approximation to the exact level density ([2]). 
Therefore, if the averaging is too accurate (see the pre- 
vious remark) the procedure leads to a curve which is 
very close to This curve remains characterized by 
strong oscillations which express shell effects. However, 
the aim of the Strutinsky method is precisely to remove 
these shell effects. 

For M = 0, the effective interval of averaging is gov- 
erned by a pure Gaussian (since Pq{x) = 1 in ([5])). In 
order to wipe out the shell effects from the averaging , 
the parameter 7 must be at least of the order of the mean 
spacing between the shells, (denoted here by hio) near the 
Fermi level: 

7 ^ (7) 

In this way, we obtain a "true" smooth density. 

If Eq.(l7|) is not satisfied, the level density remains char- 
acterized by oscillations (quantum effects) and are in op- 
position with the character of the semiclassical density. 

For M > 0, Fm is not a pure Gaussian anymore, and, 
this case becomes more complicated. Indeed, the width 
of the averaging function Fm is not really 7, because in 
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([5]) the Gaussian is modulated by the polynomial Pm- 
It turns out that to maintain this width of the order of 
mean spacing shell we have to enlarge the parameter 7 
with respect to the order M. Thus, the smoothing pro- 
cedure implies a "coherency" between these two param- 
eters. 

Usually, the couple (7, M) is determined by the so- 
called plateau condition [l8| from the typical ranges: 6 < 
M < 20, and, < 7 < 2l^. 

Due to the above " coherency" , it is easy to notice that 
the plateau is invariably moved towards the right-hand 
side (toward the largest 7), when M increases. 

Since in the harmonic oscillator the spacing between 
the shells is constant, we have in this simple case hco = 
Hlu. 

Fig.l displays the Strutinsky level density (calculated 
from Rel.®) for three values of the parameter 7.(M be- 
ing fixed). 

Curve ( a): Since 7 is too small compared to hu the 
curve (a) is characterized by strong oscillations (shell ef- 
fects) which arc close to Dirac functions (see Rcl.©)- 
Curve (b) : By increasing 7 one diminishes the magnitude 
of these oscillations. 

Curve (c): In the third case, 7 is of the order of hoj, the 
curve becomes smooth, and can be regarded as the mean 
behavior of the exact level density (the so-called smooth 
component contained in the exact level density 



F. Averaged particle-number, averaged energy, and 
shell correction 



(11) 

Where (70(e) is the exact level density given by ^ and 
Ao is its Fermi level (the last occupied level). Sometimes, 
Em.-yW is denoted as 

In Rel. (|lip . the shell correction should not contain any 
component of the smooth energy which is by definition 
already included in the liquid drop model. Consequently, 
the average density gM.-y (or the average energy Em,^^) 
must not contain any residual shell effect. 
Moreover, since 'g^^j ^is a least square approximation of 
go, one can write go ~ gM,-y + ^gM,i, with Aq ~ A. When 
the condition ((T]) is fulfilled, gM,i becomes smooth and 
the exact density go oscillates around gm,-^ making SgM,-y 
alternatively positive and negative. Consequently, from 
PT|) one can deduce formally that: 

A 

SEM.yW- J eSgMAe)de (12) 
—00 

Thus, the oscillations of the shell correction SEm.-tW ^ 
a function of A are due to the fluctuations of SgM_-y{e). 
Note: 

In practice, due to the finite size of the spectrum, 
for the shell correction, the following cut-off condition 
^last-level - A > 7 must bc required (see [l^) 



The averaged particle-number and the averaged energy 
are defined through the average density of states ^ by: 

A 

Nm^W = j gM,y(.^)de (8) 
—00 

A 

EM,yiX)= I egM,^(.^)de (9) 



The detailed expressions are given in Ref.[9|. 

In practice the upper bound of the integral giving the 
particle-number is deduced from the equation 



NM,y{X) = No (10) 

where Nq is the particle-number of the system. The 
quantity A is the Fermi level of the average density 3^,7 
( i.e. the Strutinsky level density) 

Finally the Strutinsky shell correction to the binding 
energy of the liquid drop model is defined as follows : 



III. THE SEMICLASSICAL LEVEL DENSITY 
A. Bohr's correspondence principle 

Although the Strutinsky level density is mathemati- 
cally well defined by (O, there is no physical basis for 
this smoothing. Consequently, it is necessary to build an 
other level density free of shell effects which would be 
justified by physical arguments. 

Shell effects are the consequence of the quantum nature 
of the level distribution The natural way to elimi- 
nate such effects would be to go over to the classical limit. 
To this end, we will apply the correspondence principle 
(Bohr 1923) which states that the behavior of quantum 
systems reduces to classical physics in " the limit of large 
quantum numbers" . 

Starting from this principle, the semiclassical level den- 
sity can be deduced by using the Eulcr-Maclaurin (EML) 
summation formula. 

In practice, this amounts to obtain an asymptotic series 
which contains only three or four terms, all the others 
being divergent. The first term of this series coincides 
with the known Thomas-Fermi approximation. 
The EML expansion used here is equivalent to the 
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usual standard semiclassical methods, e.g., the Wigner- 
Kh'kwood expansi on [isl , [l^ , [3 j or the method of the 
partition function [20|. The latter are based on asymp- 
totic series of powers of h. Indeed, in the correspondence 
principle, the limit of large quantum numbers amounts 
to taking the classical limit h ^ 0. 



Notes: 

In practical cases, A » cq, therefore the previous require- 
ment can be replaced with X/huj ^ 1. In this case, the 
Fermi level A must be measured from the bottom of the 
well 



B. "Asymptotic limit of large quantum numbers" 

In practical cases, the concept of "large quantum num- 
bers" must be precised by a more concrete definition. To 
illustrate this point, we start from the typical example of 
the three-dimensional harmonic oscillator. Such system 
is very simple, its energy levels are given by: 

Cn = {n + ?>/2)hhj, n = 0, 1, 2, 3, ..oo 

The quantum number n defines a shell, and huj repre- 
sents the gap between these shells. 

Thus, n is given by: 

n = (e„ — eo) here eo — (3/2) huj is the lowest level 
of the spectrum 

The correspondence principle states that for large val- 
ues of n = (e„ — eo) the quantum physics reduces 
to classical physics, in particular, the quantum (exact) 
level density go(A) defined by ([2]) should approach the 
semiclassical level density denoted here by 5sc(A). This 
can be written as : 

if n = [(e„ - eo) /M OO) then Qo ^ Qsc 

As already noted, the shell effects are mainly deter- 
mined by the small part of the spectrum {e„} which is 
located in the vicinity of the Fermi level A. Consequently 
for these levels, we have roughly e„ w A, and, this limit 
becomes: 

if n = [(A - eo) /fiw] oo, then go{X) gscW 

Since in general the Fermi level increases with the 
quantum numbers, the arguments presented for the har- 
monic oscillator are also valid for any other physical sys- 
tem. Therefore, we will consider the previous statement 
as general. However, we must now redefine hio as the 
mean shell-spacing in the neighborhood of the Fermi level 
A. As in Rel.([7|), we denote it by huj. 

if i^^^^oo, then g„(A) ^ 5,,(A) (13) 

This limit is of course unphysical. Therefore in prac- 
tice, we require the following qualitative criterion: 

if ^=i£»l then go{X)^gsc{X) (14) 
rkj 

We understand by (fT4|) that [(A — eo) /huJ] is suffi- 
ciently large compared to the unity so that 5o(A) can 
be considered as close as possible to gsc{X) with a "sat- 
isfactory accuracy" . 

Finally, the " asymptotic limit of large quantum num- 
bers" can be defined theoretically by or practically 
by Km). 



C. Two well-known analytical cases 

The procedure described above is applied in appendix 
B to two simple cases. 



1. The semiclassical level density of the harmonic oscillator 
The eigenenergies of the isotropic oscillator are : 

3 

En^,ny,n, = inx+ fly +71^ + -)hwo (15) 

nx,ny,nz = 0, 1, 2, oo 

The semiclassical level density of the harmonic oscilla- 
tor is a simple parabola (see appendix B): 



This result is well-known, and was established very 
early by means of the partition function |20|. o r more 
recently by the Wigner-Kirkwood expansion plj . 



2. The semiclassical level density of the infinitely deep 
cubic box 

For the case of cubic box with totally reflecting walls 
the spectrum is given by: 

En^nyn, = ("^ + + nj) Eq (17) 

Ef) = TT^fi,^/ (277100) , ao= side of the cubic box 
nx,ny,nz ^ 1,2, oo 

The semiclassical level density of the infinite cubic box 
(see appendix B) is also an "old " result [l^l, [2^ . 

1 TT /- 37r 1 3 11 

D. Semiclassical shell correction 

We define the semiclassical energy by a very similar 
relations to ([51fTT|) , replacing the Strutinsky level density 
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by the one of the semiclassical density. The correspond- 
ing Fermi level is denoted as A^c- 



N,,iXsc)^ J 9sc{e)de (19) 

— oo 

Esc{\sc)= j egsc{t)de (20) 

— OO 

NsciKc) = No (21) 



■^j No 

SEsc= / ego{e)d€~ / €gsc{e)de = £„ - EscjXsc) 
-oo -oo "=o 

(22) 

One must note that unlike (5i?A/_^(A), the quantity 6Esc 
does not depend on any free parameter. 



Now, one replaces the semiclassical density gsc{X+'-fX) 
by its (M + 2) first terms of the Taylor expansion around 
A (M must be even in Pm {X)). The last term gives an 
estimation of the remainder: 



y A/+1 

I — 1 

[jXf^d^^+^g^ 
^ {M + 2)\ d\''+^ 



{jX f dfeg,,(A) 
fc! dX' 



]Pm {X) cxp (-X'') dX 



It is easy to show that Pm {X) exp (— X^) behaves like 
a delta function with respect to any polynomial of order 
k < M. Consequently the first term gives back gsc{X), 
and the second has no contribution (by noticing that 
j^M+i odd). The remaining integral, i.e.: 

oo 

/(M)= / X^^+^Pm {X ) cxp {-X^)dX 

— oo 

is obtained from appendix A. 
Finally: 



IV. THE CONNECTION BETWEEN THE 
"STRUTINSKY'S LEVEL DENSITY" AND THE 
SEMI-CLASSICAL LEVEL DENSITY. 

A. Assumptions and quantitative approach of the 
asymptotic limit 

We know from the Bohr principle given in subsect. 
IIIIBI that in the " asymptotic limit" we must have 
the approximation go (A) ~ gsdX). Since the Strutin- 
sky averaging ^ gives an approximation to the exact 
level density ,9o(A), normally in this limit it should also 
give an approximation to the semiclassical level density 
(7sc(A). The role of the curvature correction term would 
be to improve the approximation. 

We start from the averaging substituting gsc{X) for 
(7o(A) in the "asymptotic limit" (|14|) . In the Strutinsky 
averaging the parameter 7 must be of the order of the 
mean shell spacing huj near the Fermi level (sec Eq.Q). 
Consequently, in Rel. pi)) we should replace fko with 7. 



9MnW - I 9sc{e)PM [ '-^^ ^ exp 



with ^ 1, and 7 ^ fiw 




Making X ~ —, we obtain: 



9MnW ~ / gsci\+lX)PM {X) exp {-X^) dX, M 



5^^,^(A)«g,c(A){l + i?,,.^(A)} 

R (\\ ^ -^^+^ Cm+'i 1 d^^+^ffsc(A) 
""^^ ^ ^ (M + 2)!g,,(A) dA^+2 



(25) 



, ,,,.,0 1.3.5 (M + 1) 

Cm+2 = i-ir^' , M even, 

with A — eo 3> 7 ^ ftw 



Eg. 1^5)1 is fundamental and gives the straightforward 
link between the semiclassical level density gsdX), and 
the Strutinsky level density g^^ ^(X) in the asymptotic 
limit A — eo S> 7, with the necessary condition of the 
smoothing procedure 7 ^ Htu. 

It should be noted that it is 'g^ ^(A) which is deduced 
from 5sc(A) and not the opposite. Moreover, gsc{X) does 
not depend on any free parameter. Therefore gsc{X) must 
be considered as the " true" smooth level density (the so- 
called smooth component of the quantum density (O), 
and A is its Fermi level. 

The quantity i?^,^ ,y(A) is the remainder of the averaging. 
From (j25p . it can easily be identified with the relative 
error: 



\Rm.i{x)\ 



5.c(A) -5„^^(A) 



5sc(A) 



even 



and represents the noise (for the density of states) of 
the Strutinsky method. In actual calculations, it is im- 
plicit, and thus unknown. It is contained intrinsically in 

5„,.(A). _ 

It is easy to check from ((25)) that the coefficient 
Cm +2 1 [M -1-2)! in the remainder decreases theoretically 
to zero as M increases to infinity (A, 7 being fixed), pro- 
vided that gsc{X) is sufficiently regular. This in principle 
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improves the average. Nevertheless, we have seen in sub- 
sec UlE] that large values of M involve necessary a slight 
increasing of 7 in the smoothing procedure, which in turn 
increases somewhat the remainder as it can easily be seen 
from Thus, it is not possible to reduce the remain- 
der without ending. 

In practical cases, the "optimal choice" M ^ 16 — 30 
leads to a very good precision, i.e. \Rj^j^^[X)\ < 0.01. 

From Rel. ((25l) it is clear that the dependency on the 
two parameter (M, 7) becomes more and more weak as 
the remainder decreases to zero. 

The only one special case where the remainder vanishes 
rigorously is that one where gsd^) is a pure polynomial of 
degree less or equal to M . This happens in the harmonic 
oscillator case. For this reason, the Strutinsky method 
must not be "judged" in this example when M ^ 2. 

Fundamentally, the Strutinsky level density appears in 
(j25p only as an approximation (and thus is not strictly 
equivalent as it is often claimed) to the semiclassical level 
density. Consequently, the smooth Strutinsky energy 
Eq.®, and the Strutinsky shell correction of Eq.pT|) 
must also be considered as an approximations to the 
respective semiclassical quantities given by Eq. (j20p and 
Eq.dlll). 

In fact, we shall see in the next subsection that the 
remainder tends also to zero with A like (7/A)*^''^^, and 
thus becomes negligible only in the "asymptotic limit" 
(7/ A) « 1. 



The relative error on the Strutinsky level 
density in two special cases 



B. 



It is instructive, to apply our result (|25p for the cases 
seen previously, i.e., the harmonic oscillator (fTC|) and the 
cubic box (|18p with the previous assumptions A — eq 3> 
7 ^ Sw. For these calculations we choose M ~ and 
M = 2 (M must be even). 

For the harmonic oscillator case, using 5sc(A) from 
Eq.([16l), wc find: 
For M = 



=(A) 



2X' 



with X ^ J ^ hiu 
(26) 

(Here of course, the mean spacing between the shells is 
constant, and we have simply hiu = hu). 
For M = 2 



For the infinite cubic box , using gsciX) from Eq. (jl8p 
we get for M = and M = 2: 



1--L2! 

16 A^ 



with 



A » 7 ^ ftw 
(28) 



5m=2,7(A) ~5sc(A) 



1 + i^^' 

^ 512 A^ 



with X ^ ^ huj 
(29) 

There also, for both cases, we obtain very similar rela- 
tions. 

Thus, in these four cases, the Strutinsky level density 
approaches the semiclassical level density, and the rela- 
tive error (remainder) tends to zero only in the asymp- 
totic limit 7/A <C 1. Moreover it is clear that in this limit, 
the Strutinsky densities becomes practically independent 
on the smoothing parameters (7, M). 

In realistic cases, the Fermi level A is fixed for a given 
nucleus. It turns out that for medium and heavy nuclei, 
the quantity A lies several units of hui above the bottom of 
the well , since hu) j, the quotient 7/ A is thus small and 
the remainder is practically negligible. Consequently, for 
these cases, the Strutinsky density of states is very close 
to the semiclassical level density. 

The relative error of the Strutinsky density for the cu- 
bic box is illustrated in Fig. 2. 

In Fig. 2a we compare the "numerical" Strutinsky level 
density (curve a) calculated by means of Eq.© to the 
semiclassical density (curve b) given by Eq. lfTS]) . Apart 
from the region near zero (very small A), and in spite of 
the "low" order M = of the curvature correction. It is 
clear, that the two densities are practically indistinguish- 
able. The theoretical link between both densities is given 
by Eq.® 

In Fig. 2b, we can see that the "numerical" relative er- 
ror {gsc — gfj^f^^-y)/gsc of the Strutinsky density (denoted 
by (b-a) /b) with respect to the semiclassical density (re- 
garded as the true smooth density), is very small, espe- 
cially in the asymptotic limit ((7/A) <C 1). In the latter, 
this error becomes close to the theoretical value 7^/16A^ 
given by Eq. 



ffM=2,7W^5«c(A)[l + 0], with A»7^fii^ (27) 

Since the semiclassical level density is a parabola, 
the derivative which appears in the remainder in 
Eg. ips)) cancels for M > 2, therefore one obtains the ex- 
act result (the remainder is 0). One should not be too 
"impressed" by this case ( see subsection IIV A[> . 



C. The relative error on the Strutinsky energy in 
the two previous cases 

The average energy Em, 7(A) can be deduced by com- 
bining Eq.® and Eg. ips]) with the assumptions of the 
asymptotic limit and the necessary condition of smooth- 
ing made in Eg. ipS)) : 
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(A) 



7 



EsciX) 
M+2 f d''+'g,,ie) 



(Af + 2)! J ^ de*^+2 

— oo 

with A — eg 3> 7 ^ ftw 



de 



(30) 
(31) 

(32) 
(33) 



where p^j,j (A) and S,^,^^ (A) are respectively the rela- 
tive and the absolute errors on the smooth (Strutinsky) 
energy £;m, 7(A). _ 

We know from Eq.f^H) that Ej^j_y{\) must be con- 
sidered as an approximation to i?sc(A). Besides, unlike 
Em. 7(A), the quantity Esc{X) does not depend on any 
unphysical parameter. As before the remainder p^^ ^ (A) 
of the Strutinsky energy must be related to the relative 
error: 



Esc (A) — Em, 7(A) 



Esc (A) 



(34) 



Once again, in order to illustrate some features of the 
Strutinsky method we apply the above result for the har- 
monic oscillator and for the cubic box with M = Q and 
M = 2. 

For the harmonic oscillator: 



Em=o, 7(A) ~ Esc{X) 



1 + ^ 
A^ 



with A 3> 7 ^ hu 
(35) 



Em=2 7(A) w Esc{X) [l-f- 0] , with X:^-f^huj 

(36) 

Where Esc{X) of the harmonic oscillator is given in ap- 
pendix B. 

In the same way, we obtain for the cubic box: 



EM=o,y{X)^ EsciX) 



Em=2. 7(A) « Esc{X) 



1 _ Ar 

16 



1 - 



512 



with X ^ J ^ hu! 

(37) 

with A 3> 7 ^ huj 



(38) 

Where Esc{X) of the cubic box is also given in appendix 
B. 

Thus, in the asymptotic limit (7/ A) ^ 1 (i.e. for 
medium and heavy nuclei), as for the density of states 
(p5)) . the relative error is small, and we have also 
Em, 7(A) ~ Esc{X). Thus, for the smooth energy, the 
Strutinsky method is a good approximation of the semi- 
classical method. A straightforward consequence is that 



the smooth (Strutinsky) energy Em, 7 becomes practi- 
cally independent on the smoothing parameters (Af, 7) 
in this limit. 

We give in Fig. 3 an illustration of the relative error 
for the Strutinsky energy for the cubic box 

Fig. 3a: Same conclusions as Fig. 2a, i.e. Em=o. 7(A) 
and EsciX) are indistinguishable 

Fig. 3b: There also, the relative error on the Strutinsky 
energy (with respect to the semiclassical energy) tends 
toward zero in the asymptotic limit (7/ A) <C 1 and ap- 
proaches the theoretical value 572/16A^ given by Eg. ([57)1 . 



D. The new understanding of the plateau 
condition on the average (Strutinsky) energy 

As seen before, the relative error plays a major role in 
the Strutinsky energy. 

From Rel.llM]), it is clear that if |Pjjj.7(A)| < 1, the 
relative variations of the Strutinsky energy Em, 7(A) are 
small compared to Em, 7(A) itself (or to EsdX)). 

For instance, if one plots Em, 7(A) as a function of 
the parameter 7 (M and A being constant) it appears 
a "plateau" in the graph. This means that Em, 7(A) is 
"practically constant at the scale of its own value" (at 
least in the interval A 3> 7 ^ ^), i.e.: 



AgA/,7(A) 

Em, 7(A) 



< 1 



(39) 



A>7^fii/j 



(which is close to ([M|)). This does not necessarily mean 
that the derivative of E^j^^iX) cancels as in the tradi- 
tional plateau condition 18']. The same remark holds for 
the Strutinsky density. 

We must point out that the relative error is propor- 
tional to the quantity (7/A)^^^^, so that the plateau is 
improved at large values of M. 

Fig. 4 shows an illustration of the plateau (defined by 
([55)1 ) for the energy in the cubic box case 

Fig. 4a: The Strutinsky energy Em, 7 of the cubic box 
is plotted as function of 7 for four values of the order 
M. The particle- number is fixed arbitrarily at Nq = 200 
with a Fermi level A = 64.255ii^o- It is clear that the 
fluctuations AEm, 7 are small compared to Em, 7- At 
this "scale" a clear plateau is noticed. 
Fig.4-b: The same figure as Fig.4-a at a "reduced scale" 
shows the important variations of the plateau, especially 
on the r.h.s of the figure. If we continue to "zoom in" 
the curve, it appears several minima and maxima (i.e. 
QEm, 7/^7 = 0). Some of them, have nothing to do with 
a plateau. For example, in the vicinity of 7 ~ 154 Eq a 
minimum occurs for the order M = 16 which does not 
really belong to any plateau. Thus our "macroscopic" 
definition of the plateau seems more " adapted" than the 
old version based on the " stationarity" of Em, 7 with re- 
spect to 7. In fact, it contains implicitly the concept of 
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the relative error which plays a central role in the numer- 
ical applications in the method. Indeed the minimization 
of the relative error (see below) avoids the ambiguity of 
the (old) plateau condition, because a stationary point is 
not necessary a plateau, whereas the minimization of the 
relative error leads indisputably to the true value (at the 
very least to the optimal value) of the smooth energy, an 
thence of the shell correction (see below). 



E. The Strutinsky shell correction 

1. The critical point of the Strutinsky method 

First, one must recall that the Strutinsky shell correc- 
tion and the semiclassical shell correction are respectively 
defined through Rel.((Tll) _and i.e. 

occupied 
occupied 

Subtracting the second equation from the first, and us- 
ing Rel. (|5D|) . one obtains a straightforward relation be- 
tween these two quantities.: 



SEM^y{X)~SE,,{X)-S^,,^ (A) 



(40) 



Since Em,^^ is considered as an approximation to 
Esc{X) (see previous subsection), the Strutinsky shell cor- 
rection SEm. 7(A) must be also regarded as an approxi- 
mation to the semiclassical shell correction SEsc{X). In 
this respect, S^^^^j (A) represents the same absolute error 
in both formulae (|30l40p . However, the essential point is 
that, this error has not the same importance in these two 
results 

Indeed, the two shell corrections SEj^j 7(A), and SEsdX) 
are obtained as the difference between two close large 
numbers (i.e. the sum of single-particle energies and their 
averages) therefore they are significantly smaller com- 
pared to these quantities (i.e. Em, 7 or Esc)- 
For example, the Ref [l^] gives for the case of ^^"^Sn 
(neutrons) the typical "realistic" values — —1122.5 
MeV, ^A/, 7 = -1132.0 MeV (the order M is not 
precised in that work), and hence, the shell correction 
^Em, 7 = 9.5 MeV is thus, much smaller than Em, ^. 
Consequently, the same absolute error Sj^j^j (A) which is 
relatively small in Eq. (|30p might become non negligible 
in Eq. ((40|) for the Strutinsky shell correction. 
In addition, in a number of cases the shell correction 
might also become so small that the relative error {in 
the shell correction) has no longer a sense. Thus, for the 
shell correction, the relative error does not play the same 
leading role as for the Strutinsky energy (or the Struti- 
nsky density), so that the (Strutinsky) shell correction 
might become strongly dependent to the choice of the 
parameter 7. This means that the error could exceed the 
shell correction itself if this error is not optimized (i.e. 
minimized) . 



The optimization of the method with the "rule of the 
relative remainder R.R.R". 



The shell correction is defined as the difference between 
two quantities ^ e,;, and Em, 7(A). Only the latter de- 
pends on the parameter 7 (and also M) through the re- 
mainder p^^ 7 (A) from (|30|) . By minimizing this remain- 
der (i.e. the relative error) with respect to 7, we make 
Em, 7(A) as close as possible to Egc, therefore we make 
SEm. 7(A) as close as possible to SEsc{X) (i.e. the true 
shell correction). Thus, the minimization of the relative 
error made on Em, 7(A), should lead to the independency 
( or at least to a weak dependence) of SEm, 7(A) on the 
parameters (7,M). Thence, we can affirm that it is the 
minimization of this relative error which is the source of 
the plateau, not the opposite. This should be the most 
appropriate way for finding the true (or the best) value 
of the shell correction. Fig. 5 displays a practical illustra- 
tion of this minimization (the so-called "rule of relative 
remainder" ) . 

We show in Fig. 5a how to find the optimal value for 
the parameter 7 in the cubic box case. We again consider 
in Fig. 5a the case given in Fig. 4a 
The steps are the following: 

• First we calculate the Strutinsky energy at the 
Fermi level Em, 7(A) as a function of the param- 
eter 7 (A being fixed), for some close values of the 
order AT. 

• We must minimize the remainder of Rel. (|34p as fol- 
lows: 



d_ 

97 



d Esc{x)-Em,^{x) 

di Esc (A) 



OEm, 7(A) 




9Sa/.7(A) 








Esc (A) 




Em,',{X) 



• Then, this quantity is plotted as a function of 7 for 
each value of M . Thence, we should look for the 
minimum of this function (relative error). 

However for a fixed A/, this function has an oscillatory 
behavior (around zero) which leads to several local 
minima. Nevertheless, due to the fact that these curves 
do not cancel simultaneously, it is possible to remove 
these unpleasant oscillations by considering the mean 
relative error over some (relative close) values of M . 
The mean relative error (over M = 16, 20, 24, 28) on the 
Strutinsky energy is plotted against 7. The minimum 
(optimal) value is found to be about ^^pt ~ 28.6£'o and 
corresponds effectively to the best value (see Fig. 5b) 
In Fig. 5b, we note a good agreement between the 
true value (semiclassical) given by the upper straight 
line Esc ~ 8093.97i?o and the Strutinsky values. 
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In fact, for the optimized 7, the calculated values 
are Ei^^ ^^^^ « 8094.52£;o,~E2o, ^^^^ « 8094.64^0, 

:B24, 7„p, 8095.21^0, ^28, 7„p, ~ 8096.31SO, for 
M = 16,20,24,28. 

At last the sum of single-particle energies is 

J2occ.states ~ 7842£:o. 

The true shell correction SE^c ~ — 251.97i?o, 
and the corresponding Strutinsky shell corrections are 
thus: ^£^16, 7„p, ~ 252.52^0, 7^^, « 252.64£;o, 

(5^24, 7„p, ~ 253.21i;o, 5E2s^ 7^^, « 254.31i;o, which are 
in good agreement with the true (exact) value. Without 
optimization the results of the Strutinsky shell correction 
will certainly be random. 



2. The semiclassical level density 

It is well-known that the level density of finite po- 
tential becomes singular at the top of the well (2^.|8|. 
For this reason, it is not possible to find a local polyno- 
mial approximation to the semiclassical level density in 
the vicinity of this singularity. In other words, the least 
square averaging ^ does not hold in this case, precisely 
because gsc{X) ^ 00 , as A ^ 0, i.e. at the zero-energy. 
This explains why the two methods do not lead to the 
same results, especially for the weakly bound nuclei. Far 
from the zero-energy there is no problem. 



3. Comparison between the two level densities 



F. Case of realistic wells: 

1. The Strutinsky level density 

First, we must note that the spectrum of the finite 
potential is composed of discrete negative levels plus a 
continuum. For this potential, the definition of the exact 
level density ([2]) must be modified by adding an appro- 
priate continuous expression gd^)'. 



It would be interesting to determine the limit where the 
Strutinsky level density deviates significantly from the 
semiclassical density. To this end, we will be comparing 
numerically the Strutinsky level density to the semiclas- 
sical ( Wigncr-Kirkwood) level density by employing the 
result of Ref [H: 



ffsc(A) = 



dNsc{\) 

dXsr 



(42) 



with: 



?o(e) 



5{e 



(41) 



For spherical potentials, the continuum is defined by 
the scattering phase shift, whereas for the deformed 
case it can be solved by the more complicated S-matrix 
method (see [1]). 

Next, one recalls that the result (|25p (which is the ba- 
sis of the present work) comes from Eq.®. The latter 
is valid for any smooth potential regardless whether it is 
infinite or not. Indeed, it is to be noted that the interval 
of averaging in ^ goes from — cx3 up to +00 so that the 
preceding demonstration remains still valid for a finite 
well. One has simply to add the continuum of (|¥T|) to 
the discrete spectrum into this integral. Thence, as for 
infinite potentials, it is clear that the Strutinsky's level 
density should also be an approximation to the semiclas- 
sical level density for finite wells. 

In practice, the rigorous treatment of the continuum 
is not easy. The standard recipe consists of using the 
discrete positive energies to "simulate" this continuum 

. These energies are usually obtained by diagonalizing 
the Hamiltonian matrix in a truncated harmonic oscilla- 
tor basis. In fact, this delicate problem seems now to be 
solved by the so-called GFOE method (Green's function 
oscillator expansion) [l^] which improves the standard 
method.. 



1 /^2to\ 



3/2 



d^j 



(A.c-f/)^/^+(f)" n 

(43) 



3,2/-, rni/2_lv2[/(A,,-C/)-i/2 



16 



This formula contains the " classical" Thomas- Fermi term 
plus an Wigncr-Kirkwood correction. In this formula 
U{~r^) is the central mean potential which contains also 
the Coulomb interaction for the protons, and fsoi/r) is 
the spin-orbit interaction. The classical turning points 
Tsc are defined by U (r^) = Xsc where Xsc is the Fermi 
level 

The numerical integration giving Nsc is made with the 
help of an "improved" Gauss-Legendre quadrature for- 
mula. 

The eigenvalues used in the Strutinsky level density are 
calculated by the code published in [20| . 
In the two methods we employ strictly the same Hamil- 
tonian and the same set of parameters, i.e. we use 
the Woods-Saxon potential with spin-orbit term and the 
Coulomb potential for the protons. For this test we work 
with the ■^''^Pb (neutrons) with a spherical nuclear shape. 
The parameters ( in MeV, fm) are: Vq = -47.083MeT^, 
Qy = 0.66Mey, Rv = 7.36/m, x = U.OMeVfm^, 
ttso = 0.55/to, Rso = 6.698/to. Their definition is given 
in M. 
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Thus, in Fig. 6, we have drawn the semiclassical level 
density (denoted by semicl.) and the Strutinsky density 
as function of the Fermi level for three cases (the numer- 
ical values of the smoothing parameters are given in the 
figure) . 

• We can check in the three cases that the Strutin- 
sky level density is practically equal to semiclassical 
density in the "intermediate" region (between the 
top and the bottom of the well) irrespectively of 
the order M of the smoothing procedure. Indeed, 
we have shown in the subsection IIVBI that in the 
asymptotic limit (i.e. for medium and heavy nuclei) 
the Strutinsky density of states (not the Strutinsky 
shell correction !) should not be very sensitive to 
both free parameters of the method. However, al- 
though with M = one obtains a good relative 
error on the density of states, one improves more 
this error at sufficient large values of M. 

• As expected, due the reason invoked above. Fig. 6 
shows that the Strutinsky densities differ from the 
semiclassical result essentially at the top of the well. 
In the "intermediate region" there is no difference. 
Near the singularity, it is more advantageous to 
choose high values for the order M . With M = 30, 
the Strutinsky calculation seem to work reasonably 
well up to about —2.5 MeV, beyond this limit the 
precision is lost. 

• However, we should not forget that it is the semi- 
classical density which is the "true" quantity. Due 
to the importance of the difference between the 
two methods near the zero-energy, the semiclassi- 
cal method must in principle be preferred for the 
weakly bound nuclei. 



V. CONCLUSION: 

Although this paper explains a number of aspects and 
"subtleties" of the Strutinsky method, we will insist on 
some essential points: 

1) The Strutinsky level density appears in the funda- 
mental Rel. (j25p only as an approximation ( and is not 
strictly equivalent as it is often claimed) to the semi- 
classical level density. Consequently, the shell correc- 
tion calculated by the Strutinsky's method should also 
be considered as an approximation to the semiclassical 
shell correction. 

2) Semiclassical quantities such as the level density, 
the energy, the shell correction must be considered as 
the "true" quantities compared to those obtained with 
the Strutinsky method. Moreover, they do not depend 
on any free parameter. 

Unlike the semiclassical method, the Strutinsky 
method contains an intrinsic noise (remainder). The am- 
biguity of the method comes from the dependence on the 
two free parameters through this remainder. 



3) It turns out that the remainder is proportional to 
(7/ A) and is defined in this paper as the relative error. 
In the asymptotic limit (7/ A) ^ 1, i.e. for medium and 
heavy nuclei, the relative error is small, especially for 
higher M. Therefore, it is found that the Strutinsky 
method gives good results for the average level density, 
and the average energy. In these cases the dependency 
on the free parameters is weak. 

On the contrary, in the shell correction the relative er- 
ror is no longer small. The shell correction might become 
strongly (7, M)-dependent. The choice of these two free 
parameters must be treated with care. In this context, 
in order to minimize the relative error, we propose the 
"rule of the relative remainder (RRR)". 

4) For realistic potentials, the semiclassical level den- 
sity admits a singularity at the top of the well. Since 
the Strutinsky method is a least-squares approximation 
to the semiclassical level density (demonstrated in this 
paper), the averaging fails near this singularity. In this 
case, the two densities are different, and it is not surpris- 
ing to note a strong dependence on the parameters (7; M) 
in this region, even if the continuum is treated properly. 
Consequently, for the weakly bound nuclei (drip-line) it 
is better to use the semiclassical method. 

Our personal conclusion, is that the semiclassical 
method with a "good numerical treatment" should, in 
theory, be quite superior to the Strutinsky method. The 
latter can be considered only as a good palliative method. 



APPENDIX A: CALCULATION OF THE 
INTEGRAL I{M) 



I{M) = 



j,A/+2p^^ (x) exp {—x^)dx, (see sub- 



secHVH) 

First, one must note that Pm {x) of (O can be ex- 

A \ gM(0) Hm+i{x) 

pressed as Pm (x) = ^mTT]^ ' ' 

Indeed, with the help of the Christoffel-Darboux for- 
mula (chap 22 of [13): 

^ Hk{x)Hk{y) Hn+i{x)Hn{y) - i/„+i(2/)i?„(x) 



E 



k=0 

one finds for our case: 



2"+in!(x- - y) 



Pm (x) = y "'"^ ' h^ (x) 



^ ffM+i {x) Hm{0) - Hm+1 (0) Hm{x) 

2^^+iM!0F(a; - 0) 
Since M is even, Hm+i (0) = 0, this gives the cited 
expression 

To calculate the integral I{M), we have to replace 

Pm (x) by the preceding result. We obtain: 

00 

Hm{0) 



I(M) 



x^'^+^Hm+i {x) exp(-a;2)dx 
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Now we use the following property (l7j ; 

oo 

J t'^Hk (st) cxp{-e)dt = 0Ffc!Pfc(s) 

— oo 

where Pk{s) is a Lcgcndre polynomial. For our pur- 
pose, we choose s ~ 1, with Pfc(l) = 1. Making k = M+1 
in the above result, one finds: 



(x) cxp{~x^)dx = 0F(M + 1)! 
so that: 

Where Hm{0) is given in the subsect Hl Dl Finally, the 
result can be cast under the following form: 



APPENDIX B: TWO APPLICATIONS OF THE 
EULER MACLAURIN FORMULA (EML) 

In the present work, to obtain some analytical results 
we employ the Euler-MacLaurin formula (EML) (l7j . 



n-l 

E 

k=l 



F{k) = J F{k)dk - \ [F(0) - F{n)] 



+ ^ [F'(n) - F'(0)] - ^ [F'"(n) - F"'(0)] + .... 

Of course this formula can be used to calculate dis- 
crete finite sums. But here, the interest of this formula 
is its application to the determination of the asymptotic 
behavior (n — > oo) of the discrete sum. In explicit terms, 
if we take a few terms (integral plus a few derivatives) 
in this formula, we obtain a quantity which is equivalent 
to the discrete sum. This means that the error (differ- 
ence between the discrete sum and its equivalent from the 
EML formula) tends to zero as n increases to infinity. 

In general, the higher orders of this formula are di- 
vergent and must be simply ignored. For this reason, 
an asymptotic expansion does not exceed three or four 
terms. 

We apply this formula for two cases: 

• The harmonic oscillator where we take just the in- 
tegral in the EML formula, neglecting thus the di- 
vergent terms (Dirac delta functions) in the result: 



D{n) 6{e-{n + 3/2)hij) 



D{n) 5 {n + 3/2) dn = gsc{e) 

inio i 



D{n) = {n+ l){n + 2)/2 is the degeneracy of the level 
n, we find: 



/ e \ ' 



1 



\huj) 4 



Now, it is easy to deduce the semiclassical energy: 
i?,e(A) « / eg,,{e)de 



• In the cubic box the degeneracy is unknown, and we 
have to evaluate a threefold integral. We work with 
a "basic" EML formula, i.c without derivatives: 



oo 

oc „ 

Y,F{k)^ / F{k)dk~^F{0) 
k=i 



Starting from that, we apply this formula in the case 
of the threefold sum: 

While all axes are equivalent, we obtain 



oo oo oo 



E E E F{nx,ny,n^) = 

F{n,j:,ny, nz)dnxdnydnz 



ria;— riy—O riz—O 

+00 + 00 +OC 





+ 00 -(-oo 



(3/2) / / F{nx,ny,0)dnxdny 





+ 00 

+ (3/4) J F(n„0,0)dn,-(1/8)F(0,0,0) 



^ h + h + h + lo (respectively) 

Here F{nx, ny,nz) = 5 (e — (n^ + n-l + nl)Ea)- 

Using the spherical coordinates we find 

' 4^3/2' ^ 8E^' ' /7' 



1 



8-Bo 



Here also we omit /o( delta function) 

3 



Therefore: 



A 



4=A=/= 







4E, 



1/2 
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